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Abstract 

We study the dependence of entropy [per lattice site] of six-vertex model 
on boundary conditions. We start with lattices of finite size and then proceed 
to thermodynamic limit. We argue that the six-vertex model with periodic, 
anti-periodic and mixed boundary conditions produce the same free-energy 
in the thermodynamic limit. We have found fixed boundary condifions such 
fhaf fhe enfropy varies confinously from zero fo ifs value for periodic bound¬ 
ary condition. We have also shown thaf fhe physical quantifies of fhe six- 
vertex model af fhe isofropic poinf does nol change in fhe case of singular 
toroidal boundary. 
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1 Introduction 


The six-vertex model has been extensively studied over the years (HI3. It was 
firstly solved under the assumption of periodic boundary conditions [[31 . Later on, 
the equivalence of the six-vertex model with periodic and free boundary condi¬ 
tions was proven by means of the weak-graph expansion [ffl . 

Nevertheless, it was also noted in a that the free-energy of the six-vertex 
model cannot be independent of boundary conditions. Moreover the six-vertex 
model was again studied in the thermodynamic limit with special free boundaries 
[|5l , anti-periodic boundaries 0, domain wall boundary ||71 [8l and recently the 
case of domain wall and reflecting end boundary conditions was also considered 
[|9l . On the latter cases of fixed boundary conditions, the free-energy and therefore 
the entropy was found to differ from the case of periodic boundary condition. 

In order to further investigate the dependence of the physical quantities, e.g 
entropy, of the six-vertex model constrained by different boundary conditions, we 
started investigating the relation among free boundary, periodic, anti-periodic and 
mixture of periodic and anti-periodic. Besides that, we have also considered the 
case of fixed boundary conditions. Our goal is to argue that the entropy of all mix¬ 
ture of periodic and anti-periodic boundaries agrees with the periodic case and that 
there are different instances of fixed boundary in which the infinite temperature 
entropy might agrees and disagrees with the periodic or domain wall boundary 
conditions. 

In particular we introduce what we called Neel boundary conditions. This case 
does not seems to be exactly solvable due to the boundary conditions, however 
we obtained some results which shows that the entropy is the same as the case of 
periodic boundary condition. We have also found that this boundary condition is 
connected with some recent generalization of alternating sign matrix ifTOl . 

The outline of the article is as follows. In section 0 we describe the six-vertex 
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model and its boundaries eonditions. In seetion[3l we diseuss the ease of periodie, 
anti-periodie and mixed boundary eonditions. The ease of fixed boundary eondi¬ 
tions is treated in seetion HI In section we have also addressed to the case of 
singular boundary conditions at the isotropic point. In appendix A, we provide 
some results for the eight-vertex model for completeness. In appendix B, we dis¬ 
cuss the Bethe ansatz solution for A > 1 and in the appendix C we discuss the 
case of separated inversions. Our conclusions are given in section |71 

2 The six-vertex model 

In this section, we introduce the six-vertex model and its partition function with 
general free boundary conditions. 

The partition function of the statistical model is a sum of all configurations 

N L 

^ = EIIII‘"A- (1) 

{£> i=lj = l 

which defines a complicated combinatorial problem. The weight can as¬ 
sume the values a(A), 6(A) and c(A), which are associated to the different vertices 
configurations of the six-vertex model (see FigurelU)- 

a a b b c c 


OJl UJ2 OJ^ U4 U>5 U>6 

Figure 1: The Boltzmann weights of the six-vertex model. 


These Boltzmann weights are the matrix elements of the so called /(-matrix 
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mm, which reads 


R{X) = 


( 2 ) 


^ a(A) 0 0 0 ^ 

0 c(A) 6(A) 0 

0 6(A) c(A) 0 

y 0 0 0 a(A) j 

One usually requires that i?-matrix is a solution of the Yang-Baxter equation, 


-Ri2(A — ^)R2‘i{X)Ri2{^l) — R2‘i{^R)Rl2{X)R2‘i{\ — /i), 

which constraints the Boltzmann weights of the six-vertex such that, 


A = 


+ 62 


2a6 


(3) 


(4) 


for any value of the spectral parameter. On the one hand this impose several 
constraints on the problem, on the other hand it allows for exact results for certain 
boundary conditions. 

The associativity of Yang-Baxter equation gives rise to the Yang-Baxter alge¬ 
bra O 


R{X - fi) (r^(A) < 8 ) r^(/i)) = (r^(/i) o TaW) r{x - /i), (S) 

where 7(4(A) = Cal{X — ^l) ■ ■ ■'Cyti(A —/xi) isthemonodromymatrix, £i 2 (A) = 
-Pi 2 -Ri 2 (A), Pi 2 is the permutation operator and A denotes the auxiliary space 
along the horizontal direction. 

The Yang-Baxter algebra is invariant under transformation of the monodromy 
matrix IfT^ . such that Ta{X) ^ GaTaW provided that 

[R{X- n),Q ®Q] =0. (6) 

For general values of A (A ^ 1), this condition implies that the matrix Q is 
diagonal or anti-diagonal matrix iTT^ . which reads 


^(0) = 1 

(1 

0 1 

|. e«=| 

( 0 

J 

1 

{ 0 

ai j 

1 1 

( 02 

0 / 
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where we are going to restriet ourselves to the ease ai = 02 = 1, whieh give us 
the relevant matriees to study the free boundary ease. The ease A = 1 has SU{2) 
symmetry, which implies that Q can be any 2x2 matrix ifT^ . 

The equations dU® provide the commutativity property of the transfer matri¬ 


ces 


r«(A) = Tr^ g^;^rA{x) 


[T«(A),T«(/x)]=0, VA,/i * = 0,1. 

However, it is worth to note that does not commute with (A). 

After a convenient representation of the monodromy matrix 


( 8 ) 


A(A) B{X) 
C{X) D{X) 


Ta{X) = 

we clearly see that the transfer matrix are simply given by 


(9) 


r(o)(A) = A{X) + D{X), 
t9)(A) = 5(A) + C'(A), 


( 10 ) 

( 11 ) 


where these are the transfer matrices with periodic and anti-periodic boundary 
conditions respectively. 

The transfer matrices (A), when multiplied successively, builds up the par¬ 
tition function of a iV x L bidimensional classical vertex model with some bound¬ 
ary conditions. 

In the case of free boundary condition, we have that the partition function ([T]) 
can be written as 


•'free 


E Tiv 


N 


0g(fc)nr(''.)(A,) 

J=1 


,fc=l 


( 12 ) 


where in each external bond one could have all possible configurations of incom 

-’{i 


ing and outgoing arrows. The matrix Qy'^ is taken from ([7]) and it stands for 


5 






the periodic or anti-periodic boundary on the vertical direction V. Its product 
(S)L. eg*’ produce all possible configurations of incoming and outgoing arrows 
along the vertical direction, where 0^ = 0,1 for A; = 1,... ,L. On the other hand, 
the product of transfer matrices for = 0,1 for j = 1,... ,N produce all pos¬ 
sible products of and and therefore all possible configurations of 

incoming and outgoing arrows along the horizontal direction. 


3 The six-vertex model with mixed boundary condi¬ 
tions 


In the homogeneous case (Aj = A,/ifc = /i = 0) the expression (fT^ is simply 
given by 


■'free 


= Tr 


V 


1 1 


(A(A) + D(A) + B(A) + C(A)) 


N 


(13) 


k=i \ 1 1 / 

L \ / k J 

Naturally, the case of free boundary contains as special cases other more con¬ 
strained boundary conditions. The case of periodic boundary condition (P) along 
both directions consist of the term 0^ = = 0, Vj, k in (fT^ . This term can be 

simply written in the homogeneous case as 


Zpp = Ttv[{T^^\A))^]=Ttv[{A{X) + D{X))^], (14) 

and was studied in |[3l . On the other hand, the case of anti-periodic {A) bound¬ 
ary along the horizontal and periodic along the vertical direction {6j = 1, 0^ = 

0, Vj, k) can be written in the homogeneous case as. 

Zap = Ttv[{T^^\X))^]=Ttv[{B{X) + C{X))^]. (15) 

This case was considered in [|6l| in the anti-ferroelectric regime (A < —1). 
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Besides the immediate eases of periodie and anti-periodie boundary ZpA as 
well the anti-periodie boundary in both direetion Zaa, one has a number of addi¬ 
tional non-trivial mixed boundaries. One can organize this in matrix such that, 


= 


^ - 2 ^ 1,1 - 2 ^ 1,2 • • ■ - 2 ^ 1 , 2 ^ ^ 

■ 2 ^ 2,1 - 2 ^ 2,2 • • • - 2 ^ 2 , 2 -^' 


^ -^ 2^,1 ^2^,2 ■ ■ ■ Z 2 N 2 L j 


2N 2 ^ 

^ free EE Zj,k- (16) 

j=l k=l 


We have that Ej ^ is the partition function of lattice N x L with (j — 1)^^^ = 
{9n, • • •, 6*1)2 “ 1)10 ^ • • •, 01)2 defining local periodic or anti- 

periodie closing along the horizontal and vertical directions, respectively. It is 
to be understood in the above notation that the binary representation of j — 1 gives 
us the vector ..., 0i} and similarly for fc — 1 and {6 *Ar,..., 6*1}. More explic¬ 
itly, for the integers j and k we have that j — 1 = 6*i2° -f 6*22^ -f • —h 6*Ar2^“^ and 
k — 1 = 0i2° -f 022^ -f-h 0 l 2-^“0 hence 


Zj,k = Tr 


V 


N 


(g)^g^™)nT(^0(A) 


m=l 


n=l 


,j = l,..., 2 ^A; = l, 


(17) 


where for a given j and k on the left hand side we have a set of (f)m and on the 
right hand side. We also have that Zi 1 = Zpp, Z 2 N i = Zap, Z 12 L = ZpA and 
Z 2 N 2 L = Zaa- 

In order to illustrate, let us consider the instance of L = iV = 2. On the 
one hand, the first column {k = 1) of the matrix (fT^ . which means (0)io = 
{02 = 0, 01 = 0)2, represents the case of periodic boundary conditions along the 
vertical direction. This is because is the identity matrix ([7]). On the other 
hand, we have the following cases along the horizontal direction: 

i) periodic boundary at j = 1 which again implies = 0 for n = 1, 2 and 
therefore 

Ei,i = Try [(T(°)(A))']; 
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ii) periodic and anti-periodic at j = 2 ((l)io = {O2 = 0, 6^1 = 1)2) 


Z2,i = Ttv [tW(A)T«(A)]; 

iii) anti-periodic and periodic at j = 3 ((2) 10 = {62 = l,9i = 0 ) 2 ) 
Zs,i = Ttv [T«(A)T(o)(A)]; 


iv) finally the anti-periodic case at j = 4 ((3)io = {02 = 1, 6*1 = 1 ) 2 ) 

Z4,i = Try [(T«(A))2]. 

For more general matrix elements, one should also take in account the Qy'^ 
boundary matrices (k > 1 ) which will reverse or not the edge configurations 
along the vertical direction depending on k value. Although we have explicitly 
given ^ 2,1 and Z 31 in the above example, in this case they are actually vanishing 
due to the property of the arrow flux along the boundary. 

More specifically, we have trivial cases which results in a vanishing partition 
function. In general, the non-trivial cases should fulfill the following rule: Let 
$ = J 2 m=i bc the uumbcr of vertical anti-periodic closings and similarly 0 = 
'Y^n=i the number of horizontal anti-periodic closings, then mod [$ — 0 , 2 ] 
must be zero. This condition comes from the parity of the flux of arrow on the 
boundary. 

In addition to that, the partition functions Zj^k is related under 7 r /2 rotation of 
the lattice with which reads 


-‘j,k 


_ ryNxL _ ryLxN 

■- ^j,k - ^k,j ■ 


(18) 


All the above features can be seen at the infinite temperature point. By tuning 
the Boltzmann weights such that a = 6 = c = 1 (A = 1/2), we are dealing with 
the infinite temperature case. Therefore, the partition function is just counting the 
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number of equally likely physical states Z = Vt. We have obtained the number of 
physical states (fl) for small sizes (u| 
are given below, 


M 2 ,2 — 


18 

0 

0 

8 

0 

10 

10 

0 

0 

10 

10 

0 

8 

0 

0 

8 


\ 


, ^ 3,3 — 


N = 

L = 

5). 

The 

cases N 

= L 

= 2 

148 

0 

0 

84 

0 

84 

84 

0 

0 

94 

84 

0 

94 

0 

0 

72 

0 

84 

80 

0 

84 

0 

0 

72 

84 

0 

0 

74 

0 

72 

74 

0 

0 

94 

84 

0 

94 

0 

0 

72 

84 

0 

0 

72 

0 

76 

72 

0 

84 

0 

0 

74 

0 

72 

74 

0 

0 

72 

72 

0 

72 

0 

0 

68 


\ 


.( 19 ) 


One can see that the case of periodic boundary along both directions allows for 
the largest number of configurations. As a result, we can show that the entropy 
of six-vertex model at infinite temperature with periodic boundary agrees with the 
case of free boundary. This is done by noting that the number of configurations 
of free boundary (Vtfree), which is equal to the sum of all matrix elements of 
Mjv,l at a = 1 /2, is larger than the number of configurations of the periodic case 
Vtpp = Zi i. Besides that, we estimate an upper bound forilfree assuming that all 
the 2 ^+^“^ non-trivial boundaries (the non-vanishing matrix elements of Mn,l) 
are equal to the maximal value flpp. This implies that. 


Opp A ^free — 2 


L+N-l^ 


PP- 


( 20 ) 


In order to take the thermodynamic limit, we raise (l20l) to the power 1/{NL), 
take the logarithm and at last take the infinite size limit. This shows that the 
infinite temperature entropy of the free boundary case agrees with the entropy 
of the boundary which allows for the largest number of configurations, therefore 
periodic boundary conditions (Spbc = Sfree)- 
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It is worth to note that this relation among six-vertex model with periodic 
and free boundary conditions at arbitrary temperatures was proven for even lattice 
sizes by means of the weak-graph expansion ||4l . Although our reasoning is lim¬ 
ited to infinite temperatures at this point, there is no restriction on the lattice sizes. 
In the next section, we are going to extend this result to different temperatures. 

Besides that, the above imbalance in the number of states among the different 
mixed boundary (fT9l) is due to the ice-rule. In the case of the eight-vertex model 
(see appendix A), this imbalance does not exist due to the additional allowed 
vertex configurations. 

3.1 Homogeneous boundary conditions 

The cases Zpp, Zap-, ZpA and Zaa are build up by the product of T®(A) or 
r(^)(A) only. This implies that we have products of commuting operators, which 
by its turn are integrable transfer matrices. Therefore, the analysis of these cases 
are much simpler than the more general cases fc. 

In order to analyze all the above cases, we will exploit some discrete symme¬ 
tries of the transfer matrices (A). 

We can define the reflection operator as 11^ = 0^=1 the parity oper¬ 

ator as = (S)m=i where denotes the usual Pauli matrices. These are 
special cases of the discrete symmetries ([7]), which implies the following commu¬ 
tation relations 


[T(°)(A),n^] = [T(°)(A),n^] =0, 

(21) 

[T(i)(A),n^] = = 0, 

(22) 


(23) 


which comes from the invariance db]) and the last relation is a direct a byproduct 
of algebraic properties of the Pauli matrices. 


10 


3.1.1 Completely periodic boundary condition 


The case of periodic boundary conditions Zpp was completely solved long ago 
by standard Bethe ansatzQ. The partition function Zpp can be written as 



(24) 


where Aj°^(A) are the transfer matrix eigenvalues Q. Therefore it is immediate 
that only the largest eigenvalue AmL(A) contributes for the free-energy in the 
thermodynamic limit, such that 


max 



where (3 is the inverse of temperature. By its turn, the largest eigenvalue of (A) 
and therefore the free-energy in the thermodynamic limit is known for all values 
of A [l3l[I]|. The entropy at infinite temperature was given asll3l 



3.1.2 Anti-periodic boundary along the horizontal direction 

On the other hand, as the case Zap lacks of the arrow conservation flux from row 
to row, it was solved much later by means of the T — Q approach [O . It was shown 
in |l6l for even L and A < — 1 that the partition function Zap produce the same 
free-energy as the case of periodic boundary condition lO . 

In order to better understand this result, it is convenient to exploit the anti¬ 
commutation rule between 11^ and T‘^^^(A). In doing so, we obtain that in the 
case that |4/) is an eigenvector of with eigenvalue A(A), one also has 

that I'k) is an eigenvector with eigenvalue — A(A). In other words, all the 
eigenvalues of are at least double-degenerate in modulus and we may 
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write, 


TW(A) = 5^Af'(A) 




( 1 ) 




( 1 ) 


A) 


( 1 ) 


(A)ff 




( 1 ) 




( 1 ) 




(27) 


j=i 

where we can make any choice of the first 2^“^ states by selecting one state be¬ 
tween each pair {A^-^^(A), — Aj^^(A)}. 

Therefore the partition function results in, 


Z^P = Trv [(T«(A))'^] 


0 odd N, 

< 2^-1 

2 ^(a5 .^^(A))^ even at, 


(28) 


which again implies that only the largest eigenvalue AmL(A) contributes for the 
free-energy in the thermodynamic limit. 


Fap — — 


lim 


ta(A!,L(A)) 


N 


(29) 


f] L,Af-Aoo LN 

The fact that Fap = Fpp establishes the following relation between the max¬ 
imal transfer matrix eigenvalues. 


lim y InAW^(A) = lim ylnA|^^^(A), 

L—>-oo 1 j L—>-oo 


(30) 


at least for A < — 1 thanks to the result in [|3. 


3.1.3 Completely anti-periodic boundary condition 

Let us now address to the case of anti-periodic boundary conditions along both di¬ 
rections Zaa- In this case, due to the relation (1231) . we have to distinguish between 
odd and even values of L. 

For L odd, the anti-commutation between 11^ and 11^ also implies that the 11^ 
eigenvalues for |^) and 11^ |^) are opposite. In this case, we choose the first 
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states to have eigenvalues +1, whieh results in 


i=i 




-i)^n" 






(31) 


Taking the traee along the vertieal direetion, we obtain 


Zaa = = Try [n"(T(i)(A))'^] 


= 


for iV odd, 

i=i 


0 


for N even, 


(32) 


whieh agrees with the matrix element values Z 2 N 2 L of (fT^ . as verified for N, L < 
5, L odd. 

On the other hand, the operators IT* and 11^ do eommute for L even, whieh 
implies the eigenveetors |\[') and 11^ |\[') have simultaneously the same 11^ eigen¬ 
values. For this reason, we ehoose an ordering such that the first eigenvectors 
alternates between If* eigenvalues +1 and —1. Therefore 

2L-1 


n^(y(i)(y))Af = 

i=i 


which implies 




+ (-1)"^^ 




(i)\ /^(i) 

(33) 




0 for N odd, 

Z^^ = Try [n"(TW(A))'^] = <J 2 ^ (34) 

2 ^ (—1)-^ ^(A^-^^(A))^ for iV even. 

, i=i 


This again agrees with matrix element Z 2 N 2 L of (fT^ . The selection rule Mod [<h — 0, 2] 
0 is naturally fulfilled by the relations (l32l) and (l34l) . 

These results imply again that only the largest eigenvalue AmL(A) contributes 
for the free-energy in the thermodynamic limit. 




(35) 
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which due to (1^ agrees with the periodie ease at least for A < —1. 

3.1.4 Anti-periodic boundary along the vertical direction 

The ease Zpa is related to Zap due to the rotational symmetry (fTSl) and therefore it 
should be related to the periodie ease. In the next seetion we will further eomment 
on this ease. 

3.1.5 Largest partition function 

Based on the direet eomputation of the partition funetion and on the above re¬ 
sults, we ean extend our previous reasoning about the agreement between the 
free-energy for free and periodie boundary conditions at infinite temperature. In 
fact, the free-energy of these boundary eonditions also agree for every set of pos¬ 
itive Boltzmann weights a, b, c, which holds true for every temperature value. 
This can be seen by noting that Zpp, Zap, Zpa and Zaa are the only matrix ele¬ 
ments of (fT^ that ean beeome the dominant term as a, b, c ehange. 

In general, we have that Zpp has the largest value for A > —1. However, the 
situation ehanges when A < —1. In the anti-ferroeleetrie phase, the largest term 
changes aeeording to the parity of the lattiee size L and N (see Table IB- 



Largest eontribution for A < — 1 

L even, N even 

Zpp 

L even, N odd 

ZpA 

L odd, N even 

Zap 

L odd, N odd 

Zaa 


Table 1: Largest element of for A < —1. 

This implies that the largest term Z ma v = max [Zpp, Zpa, Zap, Zaa], whieh 
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should satisfies 


^max < Zfree < (36) 

is the one whose free-energy F^ax agrees with the free boundary case. Therefore, 
whenever Zap, Zpa or Zaa is maximal, we have that Fap = Fj^ee, Fpa = Ffree 
or Faa = Ffree- On the top of that, the equality (l30l) imply that all of them should 
be equal Fap = Fpa = Faa = Fpp = Ffree for A < —1 and provided that they 
are allowed by selection rule. Therefore, even in the case where Zpp is not the 
largest term, we still have Fpp = Ffree- 

3.2 Mixed boundary terms 

Now we want to address to the case of mixed boundary terms Zj^^. These are 
generically a mixed product of the transfer matrices (A) and (A). As these 

transfer matrix do not commute, we loose the integrability of the statistical model. 
Although we still can diagonalize the transfer matrices exactly, we have now to 
deal with the projection of one set of eigenvectors onto another. We want to argue 
whether or not these mixed boundary partition function will lead to the same free- 
energy as the periodic case. We show that even in this non-integrable case, we can 
still extract some information about the thermodynamic limit. 

3.2.1 First row 

Although the first and last rows of (fT^ seems to be of mixed type, they are actually 
a product of transfer matrix of the same kind which still preserves the integrabil¬ 
ity. However, we are going to look at the first row case Zij in order to better 
understand the structure of the remaining mixed terms. 

Therefore, we address to the case of totally periodic boundary condition on 
the horizontal direction (0 = 0) but quite general combination of periodicity and 
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anti-periodicity on the vertical direction (<1) 7 ^ 0). The equation (fTTI) simplifies to 


Zi= Tr 


V 




m=l 




L,g ’ 


(37) 


c{4^m} 


where ( 7 ^°^ stands for the g-i\i eigenvector of T^^\\) and the function — 

We can order these eigenvectors | ( 7 ^°^) and consequently the eigenvalues (A) 

by the eigenvalues of total spin-z component, since commutes with the 

total spin operator . The eigenvalues of are given by means of Bethe 

ansatzifn. 

MA) = (a(A))^n^|^+WA))-n^|^. n = 0.....L (38) 

i=l ' A \ > 

in terms of the solution of the Bethe ansatz equation 


ci( A^) 


Qj(^Xi Aj) ^(Aj Ai) 
^(A-i Aj ) Qj(^Xj A-i) 


n 


i = 1,... ,n, 


(39) 


where the sectors n = 0 and n = L correspond to the up and down ferromagnetic 
state and the intermediate value n = L/2 {L even) corresponds to the sector with 
zero total spin. 

For the case A < 1 and L is even, the largest eigenvalue AmL(A) is non¬ 
degenerate and it is found in sector n = f. On the other hand, the for the odd L 
case, the largest eigenvalue is doubly degenerate. In this case, it has one eigen¬ 
value in sector n = and other in sector n = 

Nevertheless in both cases, the eigenvectors of the largest eigenvalues are com¬ 
plicated linear combinations of all canonical basis vectors belonging to that sector. 
The coefficients of this linear combination are all positive by virtue of Perron- 
Frobenius theorem, therefore so it is the expectation value in (l37l) . As a result, the 
action of 11^=1 eigenvector always generates a vector with at least 
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one positive eomponent in the canonical basis element of that sector whenever 
Mod[<l>, 2] = 0. Hence the coefficient is an unknown function of L, but 
most importantly it is independent of N. This implies that this coefficient simply 
do not contribute to the free-energy in the thermodynamic limit, 


Fij = — — lim 


(3 L,N^oo NL 

Here we assume that 

^ L/^gma: 


In [(ASL(A))''/itL] = - - In ASL(A). (40) 

> exp(—for large L and any 5 > 0, which was 

verified for finite lattices. 

The case A > 1 is more subtle. In this case, the leading term of the expression 
(iTTl) changes with the number of inversions <I), although the largest eigenvalue of 
the transfer matrix is always in the sectors n = 0 and n = L. This is due to the 
fact that the coefficient is non-vanishing in the sectors n = 0 and n = L 

only at $ = 0, which is the periodic boundary case. In order to see that, we write 
both and in terms of eigenvectors of such that 

1 


gW = (-1)"'^ k) 

s=0,l 

This implies that 

n I 1 1^1’ (si,..., slI , 


^) = ^(it) + (-iri;))- (41) 


(42) 


m=l 


(*> 




and since we know exactly the eigenvectors in sectors n = 0, L, we can easily 


represent them in 

basis generated by si,. 



It,. 

...t) = 

..,Sl), 

(43) 

It,. 

..,t) = 

2i / \ 

(^) 

/ L \ 

7ri^s,„ J |si,...,sl), 

(44) 


and therefore we find 


— (Amax,n=o(-^)) 


N 


2L-1 


1)'^'') + terms from other sectors. 


k=l 


(45) 
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One can see that this expression is reduced to the case of periodic boundary con¬ 
dition when j = 1 {(f) k = 0). 

The other terms in (1431) occurs for j > 1 (<h > 0). According to the selection 
rule Mod[<f), 2] = 0, the coefficient vanishes for = 1. Therefore, the 

simplest non-trivial case is = 2. In this case, as the coefficient of the largest 
eigenvalues in the sectors n = 0, L still vanishes, we need to proceed to the sectors 
n = 1, L — 1, which contain the next leading eigenvalue. Once again, we know 
the eigenvectors exactly because (A) commutes with translation operator. The 
eigenvectors of sectors n = 1, L — 1 can be written as 


vi/W ,\ = y 

m,n=l / / j 




k=l 


v/I 




( 0 ) 

m^n=L—l 




E 

k=l 


v/I 


\lk) , 


Itfc) , 


m = 1 ,..., L, 


(46) 

(47) 


whose leading term eigenvector occurs for m = 1. Therefore, we have 


Zij =-A(Ainax,n=i(A))^ + negligible tcrms in the thcrmodynamic limit, (48) 

1 j 


hence 


n j (-^max,n=l ('^)) n j (-^max,n=o("^)) PpP^ 

P L —^cxD ]_j p L —>'00 ]_j 

(49) 

whenever j contains exactly two inversions. The equality from the above limits 
comes from the fact that the ratio between the two eigenvalues tends to a non-zero 
constant when L goes to infinity, and this is sufficient the give the same limit. 

Quite generally we have the following results for A > 1. The partition func¬ 
tion Zi j is zero for $ = 2fc — 1. For $ = 2k, the partition function per site 
in thermodynamic limit can be computed from the largest eigenvalue in sector k, 
such that 

Fl^j = “4 lim 4 ln(Amax,n=A:(A)). (50) 

p L->-oo L 
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This is a consequence of the zero expectation value of nL. ef for eigenvec¬ 
tors in sectors n such that n<koxn>L — k. Specifically when $ = L, for L 
even, we have 

This result seems to differ from the case of periodic boundary conditions (l25l) 
since for A > 1 the largest eigenvalue is given AmL(A) = A„=o(A) and here 
the leading contribution comes from largest eigenvalue in the sector n = L/2 
Amax,n=-|(-^)- Howcvcr, as noted before the ratio of these eigenvalues goes to a 
constant value for large L. This confirms that the FpA = Fpp. This can also be 
proved by means of the analysis of the Bethe ansatz solution for the eigenvalue 
Amax,n=-| (•^)’ dcscribcd in the appendix B. 


3.2.2 First column 


Due to the symmetry (flSl) . the first column and the first row should produce the 
same results. However, their explicit expression are quite different in view of 
representation (fTTI) . Nevertheless, it is still interesting to look at in order to 
understand another piece of the structure of the general terms 

The partition function with periodic boundary condition on the vertical direc¬ 
tion and mixed boundary conditions along the horizontal can be written as, 

r N 1 


Zj,i = Tr 


Yl 

m=l 


(52) 


where the first and last terms Zpp and Zap were already discussed in the previous 
sections. 

In view of symmetry property (fTSl) . the result (l40l) for A < 1, the result (ISTl) 
for A > 1, we can extend (1^ as 


lim — lnA^°^ 

L->-oo L max,n= 


.(A) 


lim -lnAW^(A), 

L^oo 1j 


(53) 
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for all values of A. Note that the limit here is not taken over the maximal eigen¬ 
value of T‘^°^(A), beeause for A > 1 the maximal eigenvalue of T‘^°^(A) is not in 
sector n = ^. We have restricted at this point to the case of square lattices L = N 
with even number of sites. 

It is interesting to understand what mechanism gives rise to the selection rule 
in the case (I52l) . The product of transfer matrices in (l52l) can be written as 

N 




, ^2 




“CAT+l 


(54) 


m=l 


where given j, or equivalently {9i,... ,9^}, there exists a set of non negative 
integers kj with = N so that (l54l) is true. Hence 


N 

m=l 


= EK“T‘(a«)‘“...(a™.) 


fcjv-l 


9i 


( 0 ) 


(0) 

9i 






( 0 ) 

9n-i 


(0) I (1) 

9n-i\9n 


{9N^\9f+i) {9^+1 , 


X 


(55) 


and therefore 


%.=E(Ar)'‘(A<;')‘"...(AS.. 


{93) 


(0) \ / (0) I (1) 
9n-i) \9n-i\9n 


‘"AaS)‘"(a™.)‘"A56) 


For 0 = 1, which does not satisfy the selection rule, we may take ki = k — 1, 
/c 2 = 1 and k^ = N — k. The summation over all eigenvectors of T(^)(A) can be 
written by grouping together the contribution of pairs Ag^{X) and —Ag^{X). This 
gives us 


%I= E [ASr"A«(9r 


XX) (A‘> 


sA) {XX 


(91.92> 





- 


(57) 

choosing the eigenvectors of T^^\X) to have definite quantum number of parity, 
we see that i vanishes. 
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In order to study the free-energy in the thermodynamie limit, we ehoose the 


states 


,(0) 


( 1 ) 


to have 


gy ) to have definite projeetion of total spin- 2 ; and also 
definite quantum number of inversion operator. We also restrict ourselves to the 
case of N even and L odd for simplicity, however the final results for different 
parity of N and L are similar to the one discussed below. 

Generally, the partition function (15^ can be written as 


E (AS’)‘(AS) 

(S'oddsfi'even) 




dN-l 


fciV-1 




fcjV+1 


X 


X 


( 0)1 ( 1 ) 
^1 \92 




_j_ I" ^2+091+093 


( 0 ) 

9n-i 


(0) I (1) 

9n-i\9n 


(1)1 (0) \ / (0) I (0) 

9n \9n+i) \ 9n+i\9i 


(58) 


_j_ |'_]^^^4+a93+"S5 _|_ ^_]^^^Ar+«9jV-l+"SjV+l 


(0) 


where the numbers Oj =0,1 depending on the parity of the eigenstate 

For the case 0 = 2, we have two distinct cases. The one we called consecutive 
inversions, is the case where there is a consecutive product of two (A) matrices 
in between the product of T^°^(A) matrices, which implies fci = m — 2, /c 2 = 
2, ^3 = N — m. In the other case, we could have separated inversions, where 
the two T*^^^(A) matrices representing the inversions are apart. This means that 
ki = nil — I 5 ^2 = 1, ^3 = 1112 — — I 5 ^4 = 1, = N — m 2 , 

m 2 — mi >2, and the remaining kj are all zero. 

In order to illustrate, we discuss the case of consecutive inversions and we 
leave the case of separate inversions for the Appendix C. 

The consecutive inversions in (1581) can be shortly written as, 

2 


%i = E2(A<(A"(A»A|( 

91,92 


.sf'lA' 


(59) 


where the prime means summation over half of the states chosen between every 
pair and 11^ This case is similar to 0 = 1 but for a sign change due 

to /c 2 = 2 instead of ^2 = 1- 
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Provided we have for large L 


wmax 

,n=Ll2 pmax j 

> exp(—(5L) 

V() > 0 

A < 1, 

L 

El 

/o(o) |o(i) \ 

V^/max,?! 1 i^max/ 

> exp(—5L) 

> 0 

A > 1, 


n=0 


(60) 


as verified for finite lattices. This implies that 

L 


Fj 1 = lini 

T -M y . 


L,N^<x f3LN 


In (2^ (aSU"-"(a!;l) ) = -Fpp. 

( 61 ) 


n=0 


where the condition (l60l) assures that the projection factor {|(y'mlx,n|5'mlxy does not 
affect free-energy in the thermodynamic limit. 

As a final example of the first column i, we would like to address the case 
where we have maximal alternation of and T0)(A). In this case we have 

ki = k 2 = ... kn = I and kN+i = 0. This is equivalent to have j = jo = 
I (2^ — 1) for even N. Therefore, we have 


iV 

2 


E nASL-.AiA(9S-i|952)(9i't|9S+i)(l + (-l)‘*““” 

gi,g2,";9N-l,9N "^=1 

(62) 


where = ag^. The non-zero contributions are those where = 1 — = 




= 1 — agj^_^. This is only possible when 0 = y is even, otherwise 


every term in is zero. Therefore, we have that 


Fjgj = _ lim 


1 


L,N^oo f3NL 


lnV('2(A(°) 

\ \ max,n max/ [\i/max,nli/max/ 


n=0 


N 

i2\ T 


(63) 


Now the number of times that the projection ^(7mlx,n|fi'mlxy appears is comparable 
to N. The same assumption (l60l) tell us that the corrections are negligible and that 
the above free-energy equals the periodic case i = Fpp. 


1 +“92m + l j 
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3.2.3 General term 


By collecting all the results above, we see that the partition function is the product 
of the eigenvalues A(°)(A) and A(^)(A) and the projection of their eigenvectors 
raised to certain powers. As the leading term eigenvalues agree in the 
thermodynamic limit and the projection of one eigenvectors onto another does not 
contribute in the thermodynamic limit, we assert that Aj ^ for all values of A is 
given in the thermodynamic limit by 


^j,k 


lim 

p L —>-00 


In (A 


,(o) 




= A 


pp, 


(64) 


whenever j, k satisfies the selection rule. 

Finally, this implies that the free-energy of the all mixed boundary condition 
is equal to the free-energy for periodic boundary condition. In this context, the 
chance for a dependence of the entropy of the six-vertex model on the boundary 
conditions is the case of fixed boundary conditions, which we are going to discuss 
in the next section. 
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4 The six-vertex model with fixed boundary condi¬ 


tions 

The situation changes drastically when we consider fixed boundary conditions. In 
this case, the boundary choices was proven to produce different results for the free- 
energy and entropy in the thermodynamic limit. For convenience, we consider the 
case of square lattices L = N throughout this section. 

4.1 Ferroelectric boundary condition 

The simplest case is the ferroelectric (FE) boundary condition (see figure O. This 
boundary condition results in only one allowed physical state for any finite system 
size, 

ZpE = 1, (65) 

for any values of the physical parameters. This is a direct consequence of the ice 
rule and it was already noted in[|4ll. 


/^1 /^3 /^4 



Figure 2: The partition function for iV = 4 of the six-vertex model with 
ferroelectric boundary condition (FE). 

In fact, there are four of such boundaries associated to the four possible ver- 
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tices Ui for i = 1,..., 4. One can obtain one from another by means of rotations 
of the lattice. However, for all the four case, one has that the entropy is zero 
(SpE = 0). 

4.2 Domain wall boundary condition 

The first non-trivial instance appeared in the context of scalar products of the 
Bethe states. This is the case of the partition function the with domain wall 
(DWBC) boundary condition [fT4ll . 

Z“™C({A}, {f,}) = {^1 B(A„) ■.. B(A 2 )B(Ai) Ifr) . (66) 


/^3 /^4 ^^5 



Figure 3: The partition function for N = 5 of the six-vertex model with 

domain wall boundary condition. 

The above partition function can be written in a determinant form ifTSll . which 
was useful to the understanding of the thermodynamic limit of the six-vertex 
model with DWBC. The results for the free-energy and entropy were surprisingly 
different in comparison with the periodic boundary |I71 HI. These results and its 
finite size corrections were rigorously proven ifT^ . 
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The partition function is one of the fixed boundary conditions that one 

has inside the case of boundary Zaa- Actually, there are two equivalent domain 
wall partition function inside the double anti-periodic boundary case, which are 
trivially related by rotation of the lattice. 

The entropy at infinite temperature ~ 0.261624 is smaller 

than the periodic case. 

One could rise the question if are there other boundary conditions which result 
in an entropy different from periodic case. 

4.3 DWBC descendent 

In order to address to the question about the existence of other special fixed bound¬ 
aries, we extensively investigated the number of configurations of other boundary 
conditions. Although the complete classification of the boundary in groups of sim¬ 
ilar pattern has eluded so far, we have found some interesting cases. Surprisingly, 
we found that there is a family of boundary condition which share exactly the 
same number of configuration as the domain wall boundary. This equivalence of 
the number of configurations occurs at infinite temperature. We call these bound¬ 
ary as descendent of the domain wall boundary condition (dDWBC). In Figure IH 
we give an example for = 5. 

The dDWBC that we have been able to classify and know to persist for larger 
N can be treated by integrability tools. We use the algebraic Bethe ansatz to 
derive a recurrence relation for the , which establish a relation with the 

conventional domain wall partition function Z^^^^ of a smaller lattice. This is 
done along the same lines as [|T3]| . but here we obtain a relation between different 
partition functions and The partition function (see 
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Figure 4 : The partition funetion of the six-vertex model with deseendent 

domain wall boundary eondition, where Sj =t, | or ^ and s* is its reverse. 


Figure ID) ean be written as, 

( 67 ) 

= (S3 JJ.Ar_2 S 4 I (7^(AAr))^3_^^5(AAr_i) ■ ■ ■ 5(A2)(7^(Ai))^^^^2 |si ft'Ar_2 S 2 ) , 

where s* =t, J, or and s, is its reverse. 

In order to illustrate this, we shall take Si = S2 = S4 =|, and S3 =t as in 

ZnUM, {/^}) = {^n\B{Xn) ... B{X2)D{X,) |;tAr_i). (68) 


Using the two-site model deeomposition, where the monodromy matrix is deeom- 
posed into two parts. 


TaW 


^ 7 v-i(A) Bn-i{X) 

Cn-i{X) Ujv_i(A) 


Ali(A) Si(A) 
C'i(A) Ui(A) 


( 69 ) 
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we can derive the following relation, 


N N j-1 

{/^}) = j ('jl'Af-ll -BAf-l(AA:)-DiV-l(Aj) -BAr_i(AA;) |'ft'Ar_i) , 

j=l k=j+l 


k=l 


(70) 


where 


N 


ri = a 


(Al /^l) (^(Am /^l)) 


m=2 


_ c(Ai - /ii)c(A 2 - /xi) 
a(Ai — /ii)&(A2 — /Ti) 


r,- = 


fl(Aj_i /ii)c(Aj /Ti) 


j = 3,..., AT. 


(71) 


^ c(Aj_i - /ii)&(Aj - /xi) 

From now on in this section, we drop for convenience the indices iV — 1 of the 
operators -B(A) and -D(A), which acts in space Y{m =2 ^rn- We use the Yang-Baxter 
algebra dS]) to move the -D(A) operators to the right. In order to do so, we need the 
commutativity properties of the B{X) and the relation 


j-i 


B(v) n BiyXrri) ^ ^ (^jk 7?(Am)77(Afc), 


(72) 


m=l 


k=l m—l 

m^k 


where 


c(Aj A/j) -p-|- cx(Afc Aj) 


n 


^ ^ j, 


f^jk ^ 


b{\j — Afc) 6(Afc — Aj) 
i^k 

k=j. 


(73) 


’ 

As the state |ft) is an eigenstate of operator i7(A), we are left with the follow¬ 
ing “recursion” relation. 


N 


2 a'({A}. W) = W\Ai) 


k=l 


N 


(KAfc))^ 

j=k 


• (74) 


It is possible to obtain a concise expression for the partition function in terms 
of a determinant formula. In doing so, we replace the determinant formula for the 
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domain wall partition function IfTSlI given by, 


in the relation (1741) . We finally obtain that 



<5i 

p{^i^ P 2 ) 

p(-5i,/^3) • 

p{Xi,pn) 

Z«({A},W) = 

52 

p(A 2 , P 2 ) 

P{X2,P3) ■ 

p{\2y Pn) 


6n 

to 

Pi^NyPs) ■ 

■ p{Xn,Pn) 


where 4, p(A, /i) and / 7 v({A}, {/r}) are given by 


N 


(75) 


(76) 


6k = (-l)'+Viv-i({A} \ Afc, {/i} \ fi,)b^-\Xk)J2^j/3jk, (77) 

j=k 

c(A — /i) 


p(A,/i) = 


a(A — /i)6(A — yu) ’ 


(78) 


InUX}, {/r}) - 


Af 

n 


i<j 


(^C{jCj{bij^bjj “f" CiiCjjCiijCiji^i^CiiCjjbj^jbjj^ CijCjiduQ/jj^ 

^ PiiPjj{,^ij^jibiibjj “ 1 “ CiiCjjdijdji^ pijPji{CijCjidiidjj “|~ C-aCjjh^jhji^ 


nil {(^iibii) 


N-2 


(79) 


and we have denoted dij = a(Ai — pj) and so on. Interesting enough, the repre¬ 
sentation of the partition funetion with domain wall boundary eondition given in 
the expression (1751) with (I78H79I) holds true for all physical regimes, which means 
all values of A. 

Taking the homogeneous limit and and setting all Boltzmann weights to unity 
(the iee-point), we obtain as expeeted the number of alternating sign matrices IfTTlI . 
Therefore, this boundary produee the same entropy as the domain wall boundary 
0 . 
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The reversal of eomer arrows whieh defines the partition funetion (IFTl) eor- 
responds to exehange Boltzmann weight pairs {cJijCUs} or {u 2 ,uie} for up-left 
corner, {a;i,a;6} or {a; 2 ,Ci; 5 } for down-right corner, {a; 3 ,ci; 6 } or for up¬ 

right comer and {a; 3 ,ci; 5 } or {uj 4 ^,ujq} for down-left comer. Therefore, whatever 
configuration of compatible internal arrows is chosen, the new partition function 
is readily obtained by making the appropriate exchange of the Boltzmann weights 
of the comers. 

It is worth to note that this invariance by the exchange of the comer arrows is 
not restricted to domain wall boundaries. We can also find similar invariance in 
any other fixed boundary. 

4.4 Merge of DWBC and FE boundary 

We have also discovered another class of boundary condition which are related 
with the domain wall boundary conditions. However, in this case the number of 
physical states and therefore the entropy is smaller than the domain wall boundary. 

These boundaries are actually a mixture or a merge of domain wall and ferro¬ 
electric boundary conditions. We choose an integer number n between 0 and N. 
Starting from the upper-left corner, we fill the first n boundary row and column 
edges with arrows in the same way we would fill the domain wall boundary of type 
cus, see Figure [51 The opposite edges of these are also filled with the respective 
arrows of the opposite edges of the domain wall boundary condition. So far, we 
have used the arrows configuration of a domain wall boundary with lattice size n 
to fill our boundary of lattice size N. The remaining arrows are filled in the same 
way, but using boundary arrows of the ferroelectric boundary condition of type 

CO 4- 

Considering this boundary, we have partially frozen the arrow configurations 
of the lattice in a similar way as the ferroelectric boundary. The only differ- 
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n 



Figure 5: The partition funetion for iV = 5 of the six-vertex model 

whose boundary are a mixture of the domain wall boundary eondition and ferro- 
eleetrie boundary (fDWBC). 

enee is due to the n x n sublattiee at the upper-left eomer. This implies we are 
left with a domain wall partition funetion of size n, whieh means = 

Ylij^ nxnK^i ~ F-i) ^ Therefore, we see that the entropy at infinity 

temperature is given by 

/ \ 2 

SfDWBC = ViV/ 

For a suitably ehosen sequenee n{N), one ean obtain any value of entropy S, sueh 

that SpE ^ S < Sdwbc- 

Similarly, we ean eonstruet another kind of boundary intimately related to do¬ 
main wall boundary eondition. This is obtained by filling the first n arrows using 
domain wall boundary eondition of type oj^ as deseribed before. Additionally, we 
fill the remaining arrows with the domain wall boundary eondition of type uq, as 
deseribed in the Figure 0 This boundary eondition obtained by the merge of two 
domain wall boundary has its eounterpart in the eontext of the generalization of 
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alternating sign matrices IfTOll . 


n 


fj'l fj'2 /^3 ^'■4 /^5 



Figure 6: The partition function for = 5 of the six-vertex model 

whose boundary are a merge of two domain wall boundary condition and ferro¬ 
electric boundary (f2DWBC). 


As a result we obtain a lattice with frozen configurations, except for two square 


sublattices of sizes n and iV — n. In this case, the partition function Z 


f2DWBC ■ 


N 


IS 


a product of two domain wall partition functions of size n and N — n and the 
weights a{\i — ^j) with (i, j) outside the two domain wall square sublattices. 


r^f2DWBC 




= z. 


DWBC t 


({A}?, 


N 


I N 


Therefore, we see that the entropy at infinite temperature is given by 


(81) 


Sf2DWBC — .lina (1 — 2— (1 —;^ ) ) Sdwbci 


(82) 


whose minimum value is given by the sequence n{N) = [y] 


ofmin) _ Q 

^f2DWBC ~ 


( 83 ) 
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which is the case where we have the merge of two domain wall boundary condition 
and the ferroeleetrie boundary. This again produees entropies smaller than the 
usual domain wall boundary condition. 

4.5 Neel boundary condition 

The last ease we are eonsidering is what we called Neel boundary condition or 
anti-ferroelectric boundary. This is the case where we have the alternation of 
the arrows along the boundaries, see Figure |2l We ean find this pattern of fixed 
boundary inside the ease periodic boundary condition along both direction Zpp 
for even lattiees. There is an analogue state for odd N, which is one of the states 
of the anti-periodic boundary case Z^a- 





Figure 7: The partition function Z^^ for = 4 of the six-vertex model with 
Neel boundary condition (NE). 

In contrast with the ease of ferroeleetrie boundary condition which allows for 
only one possible state, the Neel boundary is the one which allows for the largest 
number of configurations. This is due to the ice rule, which restricts the arrow re¬ 
versal along rows and columns through the vertices ccs and ojq. The reversal of ar¬ 
rows along the lines implies in an arrow reversal along the eolumn and vice-versa. 
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As a consequence of the arrow alternation in the boundary, the Neel boundary al¬ 
lows the larger number of arrow reversals along the boundary and this propagates 
to the bulk. 

We computed the number of states of all possible boundaries for lattices up to 
N = 6 and the Neel boundary has the largest number of physical states among 
all fixed boundary conditions. Naturally, one has that the number of configuration 
of the Neel boundary ^ne smaller than the periodic boundary case flpp. Nev¬ 
ertheless, we can also estimate an upper bound value for flpp assuming that all 
possible combination 2^^ of arrows on the boundary produce the same number of 
configuration ^Ine- This provides the following relation, 

fljvB ^ Vtpp < (84) 

which implies that their entropies coincide Sne = 5'ppc at infinite temperature. 

Alternatively, one can express the partition function in terms of the mon- 
odromy matrix elements, which holds for all temperature, 

^™({A}, {m}) = (t4.... tr B(A„).4(A„_i)... D(\,)A(Xj) |t4.... tW , 

(85) 

however, it should be noted that the operators A (A), D{X) and Zi)(A)A(A) are 
not orthogonal, except in the ice-point where the product iA(A)A(A) becomes 
orthogonal. Furthermore, the partition function is not a symmetric function 
of the spectral parameters due to the fact that these operators do not commute. 
Therefore, we do not have a simple determinant formula for the partition function. 
Even though the Zi)(A)A(A) still commutes with total spin- 2 ; operator, we also 
could not obtain its eigenvalues by analytical means. 

Nevertheless, we can compute the number of configurations for the Neel bound¬ 
ary for finite lattices. We have done that for lattices up to = 20. The number of 
configuration are given in the Tabled 
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N 

number of states 

1 

1 

2 

2 

3 

7 

4 

64 

5 

1322 

6 

64914 

7 

7474305 

8 

2033739170 

9 

1305583070738 

10 

1981880443295788 

11 

7111657020627320662 

12 

60382974032926242142168 

13 

1213039653244899907872180826 

14 

57687270950680153355854587442676 

15 

6494209210696211480439308528411663853 

16 

1731204438495421321106461120147832169010790 

17 

1092829001103470428650265862752651675963745966742 

18 

1633892840599915791908254127642749411000513938128114064 

19 

5785898354977820698935460290451680551971080689572072829375890 

20 

48534629904275880189653389798729712740901732087151544103619504415896 


Table 2: Number of eonfigurations for Neel boundary eondition. 
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Figure 8 : (Color online) Entropy as a funetion of the lattiee size N for the Neel 
boundary in eomparison with the entropy for periodie boundary eondition in the 
thermodynamie limit (| 2 ^ . 

While the ease of domain wall boundary is well known to have a produet 
formula for the number of states involving faetorials, it seems that there is no sueh 
formula for the ease of Neel boundary eondition. Therefore, we are not able to 
eompute the entropy exaetly due to the faet the this boundary does not seem to be 
integrable. However, one ean eompare the Neel and periodie boundary eondition 
entropies and extraet the large N behavior (see Figure [ 8 ]). The entropy seems to 
behave as Sne = Spbc{^ — where 7 ~ 2 , whieh eonfirms the reasoning in 

The first four numbers of states for the Neel boundary eondition given in the 
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Table [2] has also recently appeared in the context of generalized alternating sign 
matrices [fTOl . In fact, there is a correspondence between number of physical con¬ 
figurations of the six-vertex model and the generalized alternating sign matrices. 
The Neel boundary condition (N even) address the problem of counting the num¬ 
ber of matrices of the following type: each nonzero entry can be either -|-1 or — 1; 
the odd rows or columns, if they have any nonzero entry, then the nonzero entries 
must start with - 1-1 and end with —1 as we fill the row(column) in crescent direc¬ 
tion of indices; the even rows or columns, if they have any nonzero entry, then the 
nonzero entries must start with —1 and end with -fl. This is to be contrasted with 
the case of alternating sign matrices, where the nonzero entries should always start 
with -|-1 and also end with -|-1 ffH. 

It was shown in [[T 0 ]| the sufficient and necessary conditions for existence of 
generalized alternating sign matrices of specific types. In the context of vertex 
models, this should correspond to the determination of the non-trivial fixed bound¬ 
ary conditions of the six-vertex model. The fact that the number of states for lat¬ 
tices up to = 3 coincides with the number of alternating sign matrices raised 
the possibility for the agreement between the number of states of the DWBC and 
the Neel boundary, which was soon dismissed by looking at iV = 41101. Whereas 
in the context of vertex models discussed here, it is clear that the DWBC has fewer 
configurations than the Neel boundary condition. 

5 Merge of DWBC and free or Neel boundaries 

In order to address to cases where the entropy varies from Sdwbc to Spec, we 
propose the merge of the domain wall boundary conditions with other boundary 
conditions which provides a large number of additional configurations. This is the 
case of Neel boundary or free boundary condition. Therefore, it appears natural 
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to propose that the merge of domain wall boundary with free boundary eondition, 
e.g. along n eentral edges (see Figure |9l). This implies that the entropy may 
vary from its value for domain wall boundary Sdwbc to the free boundary value 
Sfree = SpBc- Likewise, one could merge domain wall and Neel boundary, which 
again would produce entropies in the interval Sdwbc < S < S^e = Spec- 

F-l ^-2 /^3 F-r /^5 /^6 
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E E ^ 

ri,ri=->- 
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f4 








^5 











■ S3 

S4 

S5 





Figure 9: The partition function for = 7 of the six-vertex model with 

the mixture of domain wall boundary conditions and free boundary conditions 
along n = 3 central edges, where s,, Sj =t, | and r,, f, =^, ■(—. 


Although the conception of such examples is quite natural, the precise com¬ 
putation of the entropy for these examples has eluded us so far. 
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6 The six-vertex at A = 1 with singular toroidal 


boundary condition 


In the isotropic point (A = 1) the symmetry increases, which implies that the Q 
matrix could be any 2x2 matrix 


g=( \ . ( 86 ) 

y 921 9g22 j 

This bigger symmetry at the isotropic point was used to study general toroidal 
boundary conditions [fT9ll20]| . It was noted in [|T^ that the transfer matrix becomes 
defective when the boundary matrix Q becomes singular. This implies that the 
transfer matrix has fewer than 2^ eigenvectors. Actually, it has N+1 eigenvectors 

m, 

/_fi22\ ^ / 911 \ 

|0)^"^ = (g)( (-) (8 )(y)’ n = 0,l,...,iV. (87) 

j=l V / i i=n+l ^ ' i 

and its eigenvalues are given by 


A‘"'(A) = (911 + g2ii)[o(A)]"-"[6(A)]“, (88) 


whose degeneracy is dn = (jvA)!n! • Therefore, the effect of tuning the bound¬ 
ary parameter such that Q becomes singular has dramatically changed the Hilbert 
space and the transfer matrix has the following Jordan decomposition 


T(A) = diag(Jo, Ti, • • • , Jn) 


where Jn is the dn x dn Jordan matrix given bv [|T9l . 

^ A(^)(A) 1 0 ■■■ 0 ^ 

0 A(")(A) 1 ■■■ 0 

Jn 

1 0 0 0 ■■■ AW(A) ) 
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(89) 


(90) 






We could ask if the thermodynamic properties of the partition function would 
remain unchanged even in the case of singular boundary. In order to address to 
this point, we have to build up the partition with the full matrix Q. This is obtained 
by the successive multiplication of the transfer matrix T(A) = Tr_4 [GaTaW]^ 

N 

Z = Trv[{T{\)f]=J2dn{A^^\X)f, (91) 

n=0 

= i9u+ 922f[aiXf + biXff, (92) 

where we have explicitly used the expressions (I88H89I) . Taking the thermodynamic 
limit F = —^ limAr,Ar_).oo bi Z, we obtain that the free-energy 

= max (a, b), (93) 

which agrees with the case of periodic boundary condition [HI. Therefore, even in 
the situation where the boundary matrix becomes singular the physical quantities 
remain the same. 


7 Conclusion 

In this paper we have studied the dependence of physical quantities, like free- 
energy and entropy, of the six-vertex model on boundary conditions. 

We analyzed the case of free, periodic, anti-periodic and the mixture of pe¬ 
riodic and anti-periodic boundary conditions for finite system sizes and in the 
thermodynamic limit. We have argued that all these boundary conditions produce 
the same results as the periodic boundary conditions. 

We have also addressed the case of fixed boundary conditions and a variety of 
cases arise: 

1. We found a family of boundary conditions for which the entropy coincide 
with its domain wall boundary conditions value at the ice-point. We called this 
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boundary as descendent of the domain boundary conditions, once it has the same 
number of states as the domain wall boundary conditions and is simply related. 
We have provided a determinant solution for its partition function, whose homo¬ 
geneous limit at the ice-point produce the same number of states and therefore the 
same entropy. 

2. We have also found boundary conditions which are the mixture of domain 
wall and ferroelectric boundary conditions. This is obtained by the merge of the 
domain wall and the ferroelectric boundary conditions. There is also an additional 
possibility of the merge of two domain wall boundary conditions. 

3. Besides that, we introduce what we called the Neel boundary condition. 
Due to the direct computation of the number of states, we argued that its entropy 
is the same as the periodic boundary conditions. 

4. The entropy per lattice site of the six-vertex model assumes the values 

SpE = 0 for ferroelectric boundary conditions and = for periodic 

boundary conditions. Besides, the entropy Sdwbc = | for the domain 

wall boundary condition is in between these two values. We have showed, for the 
case of the merge of domain wall and ferroelectric boundary condition, that the 
entropy vary from its values for the ferroelectric case to the domain wall boundary 
case. In addition, we propose that the merge of domain wall boundary with Neel 
or free boundary could result in entropy values in the interval between Sdwbc 
and SpBc- Nevertheless, the explicit computation of the entropy in these cases 
is still missing. However, it is clear here that the entropy varies continously from 
SpE to SpBC- 

In addition to that, we have shown that the free-energy of the six-vertex model 
at the isotropic point does not change even in the case that the transfer matrix be¬ 
comes defective. In the case of singular toroidal boundary, the transfer matrix has 
fewer eigenstates than expected, however the free-energy in the thermodynamic 
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limit remains the same. 

Although we have found additional interesting fixed boundary conditions which 
produce entropies different from the periodic case, the complete classification of 
the boundary conditions in groups of similar pattern is still missing. Another in¬ 
teresting question would be if there exist other vertex models which depend on 
the boundary conditions. We hope to address to these problems in the future. 
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Appendix A: The eight-vertex model 

Here we summarize the results for the eight-vertex model. 

The matrix elements of (fT^ for the eight-vertex model at the infinite tempera¬ 
ture point (a = b = c = d=l) result in the number of states for all mixed bound¬ 
ary conditions. These numbers of physical states (fl) for the cases N = L = 2,3 
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are given below, 


^ 2,2 


/ 



0 

0 


0 

1 

1 

0 

0 

1 

1 

0 

1 

0 

0 

1 




, ^ 3,3 


/ 



0 

0 

1 
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1 

1 

0 

0 
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1 
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1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 
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0 
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0 
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1 

0 
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0 

1 

1 

0 

1 

0 

0 

1 

0 

1 

1 

0 

0 

1 

1 

0 

1 

0 

0 

1 


(A.l) 


where the non-zero elements satisfies the seleetion rule mod [<l> — 0,2]. In the 
general ease, one has that Qn = 2^^+^ for periodie, anti-periodie and mixed 
boundary eonditions. The eases of fixed boundary conditions, e.g the Neel bound¬ 
ary and domain wall boundary condition have the same number of states = 

states. In any case, the entropy is in the thermodynamic limit 
S = Sdwbc = Sne = In 2. This confirm the independence of the eight-vertex 
model with the boundary condition [ffl. 


Appendix B: Bethe ansatz solution for the leading con¬ 
tribution 

In order to show that the largest eigenvalue of different spin- 2 ; components pro¬ 
duce the same free-energy in the thermodynamic limit, we analyze the solution 
of the Bethe ansatz equation (l39l) . We can assume a > h + cm the ferroelectric 
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phase. A possible parametrization for A > 1 is 


a(A) = ^)sinh(A + 7 ), 

6 (A) = ^sinh(A), ^ A = cosh( 7 ), (B.l) 

c(A) = ^)sinh( 7 ), 

where g is seale faetor and 7 the anisotropy parameter. 

We replaee the above Boltzmann weights in the Bethe ansatz equations (l39l) 
and its eigenvalue expression, where we also perform the shifts A A — 7/2 and 
Aj ^ Aj — 7/2 for eonvenience. Given a, 6 , c in the region a > b + c, one can 
always obtain the corresponding 7, g, A. We checked the Bethe ansatz solution 
against the direct diagonalization in sector n = L/2 in order to identify the Bethe 
root structure leading to the largest eigenvalue of this sector. We have done this 
for L up to 10. We list the Bethe ansatz roots in the Tabled where we have chosen 
a = 2.1 , 6 = 0.7, c = 0.62. This choice for the Boltzmann weights corresponds 
to 7 ~ 0.93881. Hence the above root structure resembles that of a n = f string 
centered at —i 7 r /2 and lying along an axis parallel to the real axis 11211 [22ll . In the 
present case, this comes from the fact that whenever A* in the left hand side of 
( 1 ^ has some non-zero real part, then the ratio a(Ai — 7 / 2 )/ 6 (Ai — 7 / 2 ) becomes 
either less than unity or greater than unity. Hence for large L the left hand side 
tends to be zero or infinity. In order to have the same on the right hand side of 
( 1 ^ . one should have Xj — Aj —)■ ±7 for some i. 

In fact, the solution for the Bethe equations for L = 26 has very small devia¬ 
tion from this L/2-string pattern, see Figure [TOl 
Assuming that the exact string is given by 

A, = -if + 7 - (J - 1)) , J = 1, • • •, LI2, (B.2) 
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L 

A 

Bethe roots 

2 

3.3244 

Ai = 0. - 1.57081 

4 

7.15616 

= 0.755792 - 1.57081 

A 2 = -0.755792 + 1.57081 

6 

18.9897 

Ai = 1.24738 - 1.57081 

A 2 = 0. - 1.57081 

As = -1.24738 - 1.57081 

8 

59.0169 

Ai = 1.67202 - 1.57081 

A 2 = 0.513821 - 1.57081 

As = -0.513821 - 1.57081 

A 4 = -1.67202 - 1.57081 

10 

214.268 

Ai = 2.11256 - 1.57081 

A 2 = 0.998365 - 1.57081 

As = 0. - 1.57081 

A 4 = -0.998365 - 1.57081 

As = -2.11256 - 1.57081 


Table 3: Bethe ansatz roots for a = 2.1 , 6 = 0.7, c = 0.62. 
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Figure 10: Distribution of Bethe roots leading to the largest eigenvalue in sector 
L/2 for L = 26. There is small deviation from the L/2 string pattern. 


we obtain a reasonable approximation for kn=L /2 for large L as follows, 

Amax,n=L/2(A) ^ 

+ sinh'^(A 

= sinh^(A 

+ sinh^(A 
^ sinh'^(A 

Taking the thermodynamic limit we find 

^-fiFpA ^ (Amax,n=L/2(A))^ =a(A). (B. 6 ) 

This is the same free-energy as the case of periodic boundary in the ferroelectric 
regime for a > 6 [fT]|. 


+7/2) n 

g=-L/4+l/2 

L/4-1/2 

- 7/2) n 

<j=-L/4+l/2 


cosni^A + — 7j 

cosh(A + q'j) 

cosh(A + q'7 + 7 ) 


+ 7 / 2 ) 

-7/2) 


cosh(A + q'y) 

cosh( 7 (L /4 — 1 / 2)7 “ A) 
cosh( 7 (L /4 — 1 / 2 ) + A) 
cosh( 7 (L /4 — 1 / 2 ) + 7 + A) 


(B.3) 


(B.4) 


cosh( 7 (L /4 — 1 / 2 ) — A) 

+ 7 / 2 )e'^- 2 ^ + sinh^(A - 7 / 2 )e^+ 2 \ (B.5) 
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Appendix C: Separated inversions 

In the case of separated inversions we have 


- 


E (AS’) 


A^—m2+mi — 1 




91,92>93,34 


X (sS’Ij)”) (jS’ljf) (jf IsS’) (jS’IjS”’) [1 + 


(Cl) 


therefore ag^ and ag^ must be different for nonzero contributions. Hence in the 
above summation we should restrict {giyg^} to the sequences where ag^ = 1 — 

®93- 

For L odd, sectors n and L — n have different parities. Besides, as we can read 
from relations (|2TI) and (|2^ . for each eigenvector l^i) in sector n we have another 
eigenvector l^i) in sector L — n with the same eigenvalue. Therefore, for L odd 
we obtain 


Fj 1 = lim 


In ( 4E (AL“L.„)"^'(A£L)S(9 SL.„|sA’.>| 


L,N^oo f3NL 


n=0 


— Fpp. 

(C.2) 


under assumption of (1^ . 

For L even the sectors n and L—n have the same parity. In this case, a possible 
approximation is to retain the maximum eigenvector in sector n for gi and use the 
maximum eigenvector in sectors nil and L — n i 1 for g^. Therefore we find 

1 ( ^ 

Fj 1 = — ^ lim In I 8 


L,N^oo f3NL 


Y (A® ) 

/ j V max,n/ 


N—m2-\-m\ — l 


\ max,m / 


m2—mi —1 


n=0 L \m—n\=l 


V'^^max/ I \i/max,ri. |i?max/ \i/max,m |i^max/| 


— Fpp, 


(C.3) 


where we are making use of hypothesis (l60l) and also the following: 

lim (ASL,n) = ,lim (aJ 

L^oo ’ ' 


L —^oo 


(0) M, 
^max.ndil 


(C.4) 
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